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ABSTRACT 

The general Laguerre equations provide the simplest examples 
for testing perturbation theories. 
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PERTURBATION THEORY FOR ILLUSTRATION PURPOSES 


The aim is to use simple equations with simple perturbation functions 
which allow for simple integration of both unperturbed and perturbed 
equations. The Laguerre functions are used in this manner. 

Laguerre Functions 

Let us consider a mathematical or artificial perturbation problem 
involving the general Laguerre functions. How such a perturbation could 
arise in connection with a two dimensional wave equation with a potential 

2 2 k 

energy equal - (constant) /(x +y ) . Some of the properties of Laguerre 
functions are considered in the Appendix. The wave functions for harmonic 
oscillators correspond to Laguerre functions L with m * -% for 
the symmetric solutions and m = +% for the antisymmetric. The cases 
where m is an integer provide much simpler results when used in con- 
nection with perturbation theory. 

The harmonic oscillator solutions are reduced to the Laguerre 
functions in the following manner 


H = - + (1+0.) V , H = H q + V , 

dr 

H = --^+r 2 ,V = (2CC+0L 2 ) r 2 

o , 2 

dr 

Replacing ^ H by H } and introducing the new variable x = r 
find for the symmetric, antisymmetric solutions , respectively, 


( 1 ) 


we 


Symmetric H = - x — £ -% ^ + (1+Ob) 4 J 

dx 


( 2 ) 

2 

Antisymmetric H = - x \ + ^ 3 E n “ ^)(l+flt) 

dx 

The general Laguerre example treated as a perturbation, proceeds 


as follows: 



2 


Hy = (H+V)y = Ey 


H = - x —• - (m+1) ^ + (1+fZ.) 2 J , E n = (n+ ^)(l+a.) 


H 


x 


OL 


dx 


<-*> s + f > v - t < 1+ < K 


We shall define the unperturbed functions by 

x 1+s 
2 1-s 


oo 


t(x^s) = 


. m+1 


’ Y y n (X,S ” 


n=o 


(1-S)‘ 

corresponding to the explicit expressions 

k 


y n (x) 


V x) 


y x (x) 


= e 


= e 


= e 


= e 


x n 
2 


V /n+mNXix ) 

L vn-k/ k . 1 


k=o 


x 

2 

x 

2 

x 

2 


([ nH-l-x ^ 
[("?) - 


y 2 <x) 

From (4) the normalization integral 

* ls £7 u > 


(?) * - #] 


:< 0) = (rt-Stl) 

(3) 


(4) 


etc. 


(4 a) 


(4b) 


(5) 


whence we see that functions are normalized only in the case of m = 0 
The recurrence formula 


xy n = (2n+nrfl)y n - (n+ljy^ - (iHnOy^ (6) 

is also easily established. In the general case of arbitrary m we 
shall prefer these non -normalized functions. 
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We shall write 


Hz = Ez 


( 7 ) 


for the perturbed functions and expand them in terms of the , 

o® 

z(x) - £ c n y n (x) 
n=o 

which by means of (3) and (6) gives the linear equations 

- n f - (1 + T )C n-l + [ (n+ S I i >< 1+0L > + (2rrim+l)(y-) 2 - e] c n 

- (n+l+m) (1+ y-) c^ +1 = 0 


(7a) 


( 8 ) 


The true eigenvalues of (3) can be found from the corresponding 

determinant; however, not quite easily. These eigenvalues being known, we 

may put, say, E = E into (8) to obtain the coefficients of z ^(x) . 

We might take z D (x) = y ((1+Q6)x) . For better conformity, however, 

* & nrfl 2 

we should have to add a factor (l+<*) 2 “ 

corresponds to writing 


To the order of 06 


this 


z^(x) = (1+ ^■) m+1 ^ ((l +ot )x) 


(9) 


Then putting ( Jt ) as upper index in c^ of (7a) we have 


n-A 2 

,UL ) _ nl _ 4 > 


OL \ 2k 


n 


XI 


( 1 + 


Li ^ V (t) (-l)\n4mrfk)JL_ ( _1_ 

n-ZT / , \k/ (n-Z +k) ! (n+m) ! I . OL 

2 ' k“o 2 - 


For small 06 then approximately c 
however, easily found in this way. 
We shall, therefore, write 


a> 


( 10 ) 


= 1 . The coefficients are not. 


GO 


^z^(x)t^ = (1+ ^ L ) mf V ((!+«•) x,t) 


JL =o 


(ID 
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The transition to the harmonic oscillator, whether for non-normalized 
or normalized functions are easily performed by putting m = + % for 
symmetric and antisymmetric solutions, respectively* 

For general illustration purposes, however, the simplest case of 
m » o , with normalized eigenfunctions and self-adjoint differential 
equation, should be preferred. 

Ordinary Perturbation Theory 


(H o +V)\|r = Ei|r 

<VV* = 

* = t D + ♦ 1 

E = E + E. 

o 1 


(E-E -V)\(r 



<W»o ■ 0 < 

<W*1 - (I 1 -V)* o , 

<W*2 - (E 1 -V)* 1 + E 2 * o , 

<VV*3 - < I r v )*2 + Vl + «3+o > 

<W*4 ■ <V V) +3 + E 2+2 + Ej+j + E 4 * o , 

<W*5 - Oi-m* + Ejt 3 + Ejtj + Vl + « 5 » o 


We may take 

Ik 2 ** 



(14) 
(14 a) 

(14b) 


(15) 


1 


0 


(16) 
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Hence from the first and (n+l)st equation 

E = / Vt t i<3T • (17) 

n J o Y n-l ^ 

Using various other combinations of the above equations we find 

*1 - /V dT 

*2 ' /"W 1 " 

e 3 = /V dr - h/V'*’ <18) 

E 4 - / V *l*2 dr - E l/*l+2 dT - E 2 /V dT 

E 5 ' / V *2 2dT - E lJV dT - 2E 2 /v2 dT - E 3 TV 4 * ' 

Application to Laguerre- functions m = 0 

Looking upon the balance of this theory, considering for simplicity 
only the ground state with energy 


E = %(l+«fc) 


(19) 


we shall rather have to express it in terms of 


€ - t < 1+ f> * 

with 

- i 

E = %->[l+4e’ = %+fe-€. 2 + 2€ 3 -5e 4 + 14 6 5 --*- 



n=o 


(19a) 

(19b) 


2n - 1 


( 20 ) 
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We now have 

V = € x 


( 21 ) 


and from 


xy n - <2n+l)y n - (rri-Dy^ - -y^ , 

* Q = y o * +1 = e y l * (21a) 


we obtain 


E o = % , E x = 6 , E 2 = - € 2 , E 3 = 3 € 3 - = 2€ 3 . 

(21b) 

To find E^ and we must have ^(x) • The right hand side of the 

third equation (15) is 

€ 2 [(l-x)yi - y Q ] = 6 2 [-2y 1 + 2y^j (22) 

from the recurrence formula,, and from 


(H -E )y = ny 
v o o J n y n 


we find 


*2 



This gives 


4 4 4 4 

E 4 = -8€ 4 + 2€ +G = - 5G , 

E 5 = 25 £ 5 - 56 5 - 4e 5 - 2e 5 - 146. 5 


(22a) 


(23) 


(23a) 


(23b) 
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The energy series is easily continued into 

E = %+€-fc 2 + 2€ 3 -5£ 4 + 14 £ 5 - 42£ 6 + 132£ 7 - ••• (24) 


according to the general expression (20) • Similarly the wave function 
is^ because of 


oc 1+4C - 1 = 2+4 £ -2 aJ1+46.' 

2+OC 144*^ +1 4€ 

(25) 

= +€ - h T+4&' ] = €-2£ 2 + 5€ 3 -14£ 4 +42€. 5 -... 


and from (12) , putting JL - 0 y 

*0 

* “ Y. ( 2Sr ,ny n ’ y 0 + <«- je2+5e3 - 14e4 >yi 

n=o 

(25a) 

+ (£ 2 -4€ 3 +14€ 4 )y 2 + (€ 3 -6£ 4 )y 3 + £ 4 y 4 + ••• 


From this is seen that the successive approximations of the wave 
function according to the above method are 


♦1 = €y i * 

2 

^2 = 6 t” 2y l + y 2 3 " ( 25b > 

^3 = € 3 [ 5y l “ 4y 2 + y 3] * 

♦ 4 = C 4 £ - 14y r + 14y 2 - 6y 3 + yj 


The series can be easily continued. 

Finally it may be nice to see that the series 
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* - Y, ( 2+5r )n y n 


n=o 


really gives the correct normalization integral and the correct 
~ ^ 2 


h 2 dT ■ I 


(— > 2n 

k 2+* > 


C2±otV 


(1+ f-) 2 


n=o 


(2+OL ) 2 -«. 2 1+<X 


as it should be. 

On the other hand^ we find 


qo 

f Wdr = £ ( 2 S *) 11 [(»+%)( 1+ot) + (2n+l)(f-) 2 - (2n+2)(f") : 


n=o 


d+f-) 2 

1+ OL 


OU2 


OO 


0<i +a -> - <T> 2 ] + <->Y n <2^r>" 


Z_. 

n=o 


(1+ T ) r, ,ol . 2 . ,oc.2~i , (1+< 2’ ) 


1+OC 


[%(1+0C) - (f-) 2 + (Sf-) 2 ^ = ^(1+OL)- 


l+oc 


Hence y 


(26) 

energy 

(26a) 


) 

(26b) 


E 


%(!+«-) 


(27) 
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APPENDIX 

Additional Considerations Concerning Laguerre Functions 
Consider the two-dimensional wave equation 

>2 „ \ i j p 



+ ^ + E + 

Sy 


♦ - 0 , r = J 


x 4y 


in polar coordinates 


I ai + 1 a_ + i. 

\> 2 r r 2 


sf 2 


+ E + -f t = 0 


Make the substitution 


(Al) 


(A2) 


f = y(r)e 


im^ 


(A3) 


to obtain the radial equation 

,2 


f d 2 , 1 d m 2 . „ . z\ , . „ 

^ + 737 --2+E+ 7 Jy(r) - 0 


Putting 


n 


-I 


m 

r? 


+v ~«J -E* r 

e 


V“0 


it is found on requiring the series terminating on n that 


(A4) 


(A5) 


Z = a/-E (2n+nri-l) 


(A6) 


or 


E = 


4(n+ — j - ) 


Therefore, making the transformation 


(A7) 


r = 2 'n/-e' x , x = r/2 n/-E ' 


(A8) 
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we may write 



which for integral n has the solutions 


y 


x m 

2 2 i m /• x 

x L nh» (x) 


0 


> 


the L m . being the well-known generalized Laguerre polynomials 
n+m 

degree n . 

The ordinary Laguerre polynomials L^(x) obey the equation 


^~2 + ^ 1_x ^ dx + L - = 0 J 


n 


or 


{■ 


r + ( 1-x) 1- n+m l L = 0 

, 2 dx I n+m 

dx 


hin^ 


On differentiating m times the latter equation we obtain 

,2 


I- 


dx 


— "i 

nl ] 


+ (m+l-x) — + n ^L^(x) = 0 


n+m 


m 


w*> - 0 


Substituting (A10) into (A9 ) > we have 

x 77 + (nrf ' 1 - x) k + n 

dx 

m d m 

the same equation if L , (x) = — =r L , (x) . 

n+m n+m 

The polynomial may be defined as 


t 


m 


W x > = x e 


m x 1 d n n+m -x 


n; , n 
dx 


x e 


from which their explicit expressions 

n 


^< x > - 1 (s) 


k 


k=o 


(A9) 

(A10) 

of 

(All) 

(Alla) 

(Allb) 

(Allc) 

(A12) 


(A12a) 
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as well as their differential equation (All) are easily obtained* 
Another equivalent and often useful definition is 


xt 

1-t 




(1-t) 


nrt-1 


o O 


■I 


T ® / \ 


n=Q 


(A12b) 


It should be observed, however, that by this definition 


0*> ■ 


(A12c) 


for integral m . For non-integral m , we must stick to the above 
direct definition. 

Example: 2 345 

L^(x) = 1 - 5x + 10 ~ - 10 If + 5 47 " 57 9 


L^ 3) (x) = - 10 + 5x - jj = <- 1 ) 3l 2+3 (x) 


(A13) 


If in (A9) we make only the substitution 
m 

y = x 2 Y 


(A14) 


we find the equation 


Y = 0 


f d 2 . , - , d x , m+l~ \ 

V d? ( ) dx " 4 + n + ~TJ 

referred to in the former Eqs. (2) and (3). 

If in the function Y we replace the argument x by x(l+<L) 


(A14a) 


Y((l+0Ox) = Z(x) 


(A15) 


we^ of course, obtain the equation 

f 2 

i A 2 ^ V dx ” 4 


x + (nH-l-x) — - (1+Qt) 2 jr + (n-f — jpHl+dt) V Z = 0 

dx J 


A 


(A15a) 
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The normalization factors 

eo 


/V* l* wt', (x) X m d X - te&P- 

J rrhn n + m v 7 nj 


are most easily found from the integral 


00 

J* x m e X \Kx, s)\jr(x, t)dx 
o 


rr 

(l-s) nri ' 1 (l-t) ,ttfl X 


Z (nH-m)J n n 
n! ' S 

n=o 


(Ai6> 


m.' 


(1-st) 


trrt-1 


(A16a) 



